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ABSTRACT 



A theoretical and experimental investigation of the stress dis- 
tribution and deformation of twisted blades under combined tension and 
torsion was conducted by LCDR Lawrence M. Hogan, USN, at the Aeronau- 
tical Engineering Department of the U. S. Naval Postgraduate School, 
Monterey, California. Theoretical expressions were developed to relate 
the deformation of a twisted "blade" of rectangular cross section and 
constant rate of twist to an applied load and moment. These expressions 
were found to be linear for the rates of twist investigated and depended 
only on the blade material and geometry. The theoretical strain distri- 
bution was found to be essentially parabolic across the chord, with the 
shape determined by the initial twist rate and the angular deformation. 

The experimental investigation was conducted with three specimens 
of constant cross sectional aspect ratio but different rates of twist. 
Correlation with theory was found to be excellent for the chordwise dis- 
tribution of strain and generally good for the linear relationships 
linking the deformation with the applied loads. Recommendations were 



made for areas of further investigation. 
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1 . 



Introduction* 



An analysis of the stress distribution and the distortion of twisted 
blades subjected to combined tensile and torsional loads has several fun- 
damental engineering applications. Propeller and turbine blades are two 
examples of the special case in which the blade is under tension but es- 
sentially free of torsion. High rotational velocities place these blades 
in tension due to centrifugal force while only an insignificant amount of 
torsion results from the aerodynamic loads on the blade. The need to 
predict the resultant stress distribution and accompanying distortion 
of these blades is evident. 

This paper is an attempt to evaluate the applicability of a simple 
theoretical approach to this problem. Specifically, the analysis involves 
a twisted "blade" of constant rectangular cross section and constant rate 
of twist, subjected to known applied torsional and tensile loads. The 
predicted results for three different rates of twist are compared with 
experimental results obtained from three test specimens. 

The strain at any chordwise station is determined from an analysis 
of the geometry of the twisted blade. All deformations are restricted 
to the elastic region. The central chordline of the blade passes per- 
pendicularly through the longitudinal axis. A key assumption is that 
this chordline, originally straight, remains straight in the deformed 
state. Thus the displacement of a chordline a small distance from a 
reference chordline is expressed as a linear displacement in the longi- 
tudinal direction and a rotation about the blade axis. The strain at 
any chordwise station is established as being dependent on ititial rate 
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of twist, on the longitudinal strain at the centerline, and on the 
angular deformation. 

The basis for the theoretical analysis can best be described as 
modified membrane theory. Expressions are obtained for an applied ten- 
sile force and for an applied moment in terms of membrane stress and 
the associated geometry of the blade. The stresses are then expressed 
in terms of the previously determined strain relationships by means of 
Hooke’s Law. This allows the applied force and moment to be expressed 
in a form linearly dependent on the longitudinal strain at the center- 
line, e^, and on the non-dimensional angular deformation, ")C . These 
solutions are modified by the torsional relationships for a twisted pris- 
matical bar of rectangular cross section. The resultant linear expres- 
sions are then inverted to solve for e and X in terms of the applied 

o 

tensile and torsional loads. 

Test data was obtained to compare the actual relationships with 
those predicted by theory for two special cases of applied load and 
moment: 1) An applied tensile load with zero applied moment and 2) 

an applied moment with zero applied tensile load. 
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2. Analysis. 



A diagram of the model to be analyzed is shown in Fig. 1. The 
model is considered to be an axisymmetric helical blade machined from 
a circular rod. The angle of twist, represented by 0 in view (a), varies 

dQ 

uniformly over the length, L, giving — ■= constant. The edge of the 
blade describes a helix, the projection of which is a sine wave as is 
evident in view (b). Any cross section of the blade normal to the 
z-axis is essentially a rectangle of length 2a and width r .h. 

The helix angle, 0, describes the angular difference between the 

tangent to the helix at any radius, r, and a line parallel to the z-axis, 

as shown in view (c). Thus 0 varies linearly with r from zero at r * 0 

to 0 at r =s a. 
max 

As a first approximation in the analysis, the blade thickness is 
assumed very small so that it approaches a pure membrane. The blade is 
therefore considered to have zero bending stiffness and the applied load 
must be taken by tensile stresses in the center surface of the blade. 

An inspection of the applied load, P and T, shows the stress in the 
radial direction to be zero. Therefore, the significant stress for this 
analysis is that represented by , the stress directed alopg the tan- 

gent to the helix. Considering then the element of cross section surface 
on which C T^_ acts, tdr, the elemental membrane tensile force is repre- 
sented by cr^tdr, view (c). The component of the tensile force in the 
z direction, CT^tdr cos 0, integrated over the chord must be equal to 
any applied tensile load, P. The component of the tensile force perpen- 
dicular to the z-axis, C7“ t tdr sin 0, introduces a moment about the 
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centerline since it acts at a distance r from the axis. Thus the inte- 



grated moment over the chord must be equated to the applied moment, T. 
These relationships are expressed in eqns. (1) and (2). 



The assumption of a rectangular cross section is not exact in actu- 
ality due to the machining process. The blades herein analyzed were 
machined from round stock. The stock was mounted in a manner so as to 
pass spirally by a cutting bit, as illustrated in Fig. 2. If this 
cutting bit were of infinitesimally small radius, it would act only at 
a particular z -station and an exact rectangle would be formed in the 
cross section normal to the z-axis (neglecting the arc on the ends). 

The finite radius of this bit caused it to machine on adjacent z -stations 
as well, introducing a variation in the thickness across the cross sec- 
tion. Thus the thickness t can be found from Fig. 2 to be 



where R is the radius of the cutting bit. 

Rearranging and defining the non-dimensional cutter bit radius, 
T' , as R/h gives 




a) 




sin 0 



( 2 ) 



t = h cos 0 - 2R(1 - cos 0) 



(3) 




( 4 ) 
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Fig. 3 shows an elemental length of blade, dz. Line OP is the 
helical path at radius r over the elemental length. The circumferential 
displacement over the length dz is seen to be rd0. Thus the helix angle, 
0, is defined by 



0 



tan 



-1 



rd6 

dz 



(5) 



The non-dimensional rate of twist, K, 
jf , are defined as 




and the non-dimensional radius, 

( 6 ) 

(7) 



Thus eqn. (5) becomes 
0 = tan ^ K f 



( 8 ) 



giving 



0 = tan K 

max 



sin 0 * 






COS 0 =: 



VTTl^p- 



(9) 

( 10 ) 



(U) 
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/ 



( 12 ) 



1 - cos 0 = / - 



V/ ± - / 






The radical in the numerator of eqn. (12) can be expanded by the 
binomial theorem to give 



V 7 * * 2 / 2 = / * i k- 2 ^ 2 - 7 



(13) 



Dropping higher order terms there is 



1 - cos 0 - 



jLiLeL 



V / * 



(14) 



Thus eqn, (4) becomes 



t = h 



/ - rK a f z 

Vl 



(15) 



Substituting these expressions back into eqns. (1) and (2) and 
changing the variable of integration from dr to d^ , produces 



T =2 



2 ah I CL —r de 

J o t / * ^ 

a'h y <r t 



'f 

n‘(l-r K'}*) 

/ / k 2 / 2 






(16) 



(17) 



Now an expression for (7~ t in terms of deformations will be de- 
rived, Referring once again to Fig. 3, an important simplifying assump- 
tion is made which claims that the radius SP, originally a straight line 
normal to the z-axis, remains a straight line normal to the z-axis in 
the deformed state, S , P f . 
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The unstrained length of an element along the tangent is 



OP = — - dz |//V K 2 f Z 
COS 0 J 



(19) 



The strained condition is defined by a change in length parallel 
to the z-axis, du, and an angular change, d 5^ . The length along the 
tangent to the helix in the strained condition, OP', is found as follows 



OP' 


Hi 

II 


+ du.) 2 ' 4- [ 


y(d& + dr) ] Z 


OP' 


. d* V 


(U cl') 2 - + 


[Ke'+ r')] 2 


OP ' 


1 - di \ 






Replacing u 1 and by the non-dimensional values e Q and /K/K, 


tively, 


and expanding gives 




OP 1 


' = dz V/ 


4 2 e„+ e 2 


+ K 2 f 2 ( /*- 2. X /k + 2/V) 


OP 1 


1 = d Zr y 


(UK 2 f z ) 


+ [ej2+e*) + ZKf z (2. f %0j 


OP’ 


1 - d? 




e 0 (2 + e o ) 4-ZKf z (Z + */ K ) 

/ j- Ir 2 j? 2* 



( 20 ) 



The strain along the tangent is now expressed by e^. 



OP' - OP 
e t “ ’ OP 
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( 21 ) 



e 

t 




e„ ( 2 4 e. ) + + V K ) 

/ 4 K z f* 



- / 



1/2 

The radical may be expressed as the binomial expansion of (1 + F) 
and substituted as follows: 



- ( / F -jF i 



) - I 








+ • - - 



( 22 ) 



Eqn. (22) is therefore the strain along the tangent to the helix 
at any radius, ^ . An order of magnitude analysis of the F term is 

now desired for simplifying purposes, e q will be very small compared 
to unity, as will be %/K if "large" initial twist rates are consid- 
ered; i.e., if the ratio of angular deformation to initial twist is 
small. Neglecting these terms, F in eqn. (22) becomes 



F 



1 ~ 



e 0 » x,Kf* 

/ f K 2 ^ 2 - 



(23) 



The denominator of F^ is of the order of magnitude unity while the 
numerator is of much smaller order. It is therefore possible to neglect 
powers of F^ and the expression for e fc becomes 



e 



t 




e„ +■ XKf* 



(24) 



Eqn. (24) is the simplified expression for the strain along the tangent 
at radius £ . It can be rearranged to produce a linear relationship 
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2 2 2 

between e^(l+K ^ ) and ^ of the form, y a b + mx as follows: 

e t (1 + K 2 f 2 ) a e Q + (K %) $ 2 (25) 

Since the strains are limited to the elastic region. Hooked Law 
may be applied to give 






- E e t 



cr t - E 



e a + KX f 



Z cZ. 



J 4- K z f 



(26) 



Eqn. (26) substituted into eqns. (16) and (17) produces 
P . df (27, 

J, V + K 2 -^) 2 - 

t- Jr (28) 

J 0 (/^^ Z f 2 ) z 



Non -dimens ionalizing these expressions there is 



p* * 



p* _ 



T* = 



2 G-h E 
./ 

T 



V (/+>r z f 2 ) z f (/ + 



2 of- h E 



r / — — f di + % c/f (30) 

^ + i U+K'f 1 -)* T 



t* . e. 
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Eqns. (29) and (30) are the fundamental theoretical equations to 
be solved. The various integrals are constants determined by the non- 
dimensional rate of twist, K, and are redefined as coefficients 

so that 



pvc , , Go X/ 



(31) 



T* = e o + K 



(32) 



where c< u = 

These expressions can be inverted to solve for e Q and . 



e = p* . 7-tl* r * 

° <=^ 2 2. " ^i 2- * aK 'fa?" 

H 

■y _ Z^J L l n* , 2ili 7-* 

^2 ~ °<IZ 



(33) 



(34) 



The coefficients of P* and T* in these equations are also constants and 
are redefined as @ coefficients. Hence 



e = 
0 


ft. p* * / 


In T* 


(35) 


X- 


h,P* * 1 


hz T* 


(36) 


where |3 


12 = 21* 







Solving for the o< coefficients in terms of p coefficients gives 
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fill 


11 


@11 (Izi ~ fiiz fill 


0^12 = 


~ fill 


fin fizz - fi tZ fill 


°^21 = 


- fix I 


fill fizz " fill fii.1 


<=<22 = 


/fin 


fin ftziL ~ fi 1 



(37) 



Eqns* (35) and (36) give and X for any combination of applied 
P* and T* if membrane theory alone adequately defines what is taking 
place* This solution will contain an appreciable error if the bending 
stresses, which the theory assumes to be zero, are of sufficient mag- 
nitude* Any such stresses would have a negligible component in the di- 
rection of the applied tensile load but they may constitute a significant 
portion of the torsional component. If it is assumed that bending 
stresses have a negligible effect on P*, membrane theory alone will 
adequately describe P* as in eqn. (31). If the cross -coupling coeffi- 
cients, o< ^2 an d 2 i’ retain their equality in eqns . (31) and (32), 

then the only significant departure from membrane theory would be 
evident in °^22* 

Assuming C< and ? i are a H constants as determined 

by membrane theory, a different theoretical basis is chosen for 
evaluating c>< ^ • 

From eqns* (34) and (36) there is 



p 



o<' II 



22 



2 z i : 



(38) 
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o< 22 can written as 



2 . 2 . 

_ °^n ^ 12 , /Z. 

°< 22 



22 = 



' 7 ^ + 



* /Z 



( 39 ) 



In the theory of torsion for prismatic bars of narrow rectangular 
cross section the angular deformation is related to the applied 
torque by^ 






T 
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I 

Cr h ^ ( • (p (?7 Q. - .21 h) 



(40) 



A recent investigation by Chen Chu has modified this relationship by 
a factor, f , so that it applies to rectangular bars with an initial twist , 



I 

^ fQrlr)3( t (j>(fi7cL-,2~lh) 



where 



f = I (/+/.)(’ </>„*, ‘ 



^Timoshenko, S,, and J. N, Goodier, Theory of Elasticity 
(New York: McGraw-Hill Book Company, 1951) p. 278, 

2 

Chen Chu, The Effect of Initial Twist on the Torsional Rigidity 
of Thin Prismatical Bars and Tubular Members, Proceedings of the First 
U, S, National Congress of Applied Mechanics , (New York: The American 

Society of Mechanical Engineers, 1952) p. 265. 



(41) 



(42) 
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Non-dimensionalizing eqn. (41) there is 



1L . 



2a s E 



T* 



&h l (.U7a.-.2lh)[l + (I +/*)($„*> ( ^ ) Z ] 



( 43 ) 



For the condition of zero applied tensile load in eqn, (36), 



0 ,2L 

P 22 T* 



(44) 



Eqn. (43) is therefore an expression for ft 22 ^ ase ^ on t ^ ie theory 

of torsion. Since this relationship is for a rectangular cross section, 

a slight modification of the h term is in order for the present case. 

As a first approximation it is assumed that the average blade thickness 

normal to the z-axis, h , may be substituted. 

* avg ’ J 



avg ~ J cos 0 






avg 



-/ h{l-rK % f l )df = h(l - 

''Vi 



(45) 



(46) 



Substituting into eqn. (39), ©< ^ is found to be 



cK 
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__ g 



n 

(47) 

This establishes the last unknown influence coefficient in eqns. 
(31) .and (32), allowing a complete determination of these equations as 
well as eqns. (35) and (36). These equations now completely describe 
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the deformation of a twisted blade under any combination of applied 
tensile load and torque, within the range of validity of the previously 
listed assumptions. 
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3. 



Test Equipment and Procedures 



Purpose 

The purpose of testing was to obtain experimental data for longi- 
tudinal and angular deformations for known values of applied torque and 
tensile load utilizing blades of different initial twist rates. From 
these relationships experimental values were available for: 

a. Young* s Modulus, E. 

b. The strain distribution across the chord for comparison with 
the parabolic distribution predicted by eqn. (25). 

c. ^ coefficients in eqns. (35) and (36) for comparison with 
the theoretical values. 



Background 

The basic requirement for testing a twisted blade under the con- 
ditions described in the Analysis is to have a test apparatus which 
allows the application of arbitrary, but known, external tensile loads 
and moments. Lacking any such apparatus, tests were devised to inves- 
tigate the two extreme situations: 

Test I: An applied axial load with no applied moment. 

Test II: An applied moment with no applied axial load. 

Test Specimens 

Test I can not be accomplished in a normal tensile test since any 
tendency for the blade to untwist in reaction to applied tension would 
be resisted by the gripping head, thus introducing an undesired applied 
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torque of unknown magnitude. Even a "swivel" head would tend to intro- 
duce such a torque due to the large frictional forces involved. For this 
reason the test specimen was machined to twist in one direction over 
half the length of the specimen and then the twist was reversed over 
the other half. These halves were therefore symmetrical except for the 
reversed twist. Any tendency for one half to untwist under tension was 
exactly equal to that of the other half. The center of the specimen was 
then free to untwist and, since one twisted section had no effect on the 
other, each half was essentially a twisted blade in tension but free of 
applied torque. 

The test specimens were machined with a 1.25 inch diameter cutting 
bit from one and one-half inch diamater 7075-T6 extruded aluminum rod. 

The design was such that a cross section normal to the z-axis was 
essentially a rectangle of width h = .2 inches and length 2a s 1.5 inches. 
Thus r = R/h = 3.125. 

Three different test specimens were machined, as shown in Fig. 4. 

The machined twist, as characterized by dO/dz, was 22.5, 11.25, and 
zero degrees per inch for specimens #1, #2, and #3 respectively. At the 
center of specimens #1 and #2, prior to reversing the twist, a short 
section of the rod was left intact to allow for the measurement of the 
deformation twist, V . The twisted sections were approximately nine 
inches long so that strain readings taken at the midsection would be 
free of end effect. 
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Longitudinal Deformation 



Longitudinal deformations were measured by means of strain gages, 
type FAP 12-12 (S13) and type A-7, bonded to the specimen at the center 
of each machined section. Placement was such that strains along the 
tangent to the helix were obtained at jf = 0, 1/2, 3/4, and 1 on each 
side/edge of the blade. This duplicate placement plus the “double- 
ended" feature of the two twisted specimens allowed for a good 
multiplicity of data. 

The strain gages were fed through a 20 channel Baldwin-Lima- 
Hamilton Switching and Balancing Unit and into a Budd DATRAN Digital 
Strain Indicator, shown in Fig. 5. This arrangement combined the 
advantages of a direct strain reading in micro inches per inch with a 
rapid reading capability on all gages. 

Angular Deformations 

If a circular disc of radius 11.46 inches were attached as a 
collar to the specimen, a radial angle of five degrees would subtend an 
arc length of one inch on the circumference of the disc. An angle 
indicator, shown in Fig. 6, was devised with this principle in mind. 
Rather than the complete disc, a three inch width of the disc was 
clamped around the specimen and the circumferential arc was calibrated 
in inches, A sighting gage fixed in a suitable stationary position then 
allowed circumferential linear readings for angular deflections. These 
readings of ^ , taken in sixteenths of an inch, were converted to X 
by the multiplying factor 
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lb 57. 3 




( 48 ) 



where L is the length of the twisted section. 



Test I 



Test I was conducted in a 300,000 pound Riehle Testing Machine, 

The specimen was positioned to accept a tensile load and gripped in 
fixed clamping heads at either end. The angle indicator was clamped 
between the two twisted sections for specimens #1 and #2, with the 
sighting gage appropriately placed to read instantaneous angular 
deformations. The test set-up is shown in Figs. 7 through 9. 

The strain gages and angle indicator were zeroed at a preload of 
1000 pounds, referred to as M gage zero". Readings were taken at 2000 
pound intervals to 10,000 pounds gage. This range of loads allowed 
the maximum strain to remain well below the elastic limit of the material. 



Test II was a torsion test with P s 0. For this test a torsional 
apparatus was set up as shown in Figs. 10 through 12, The twisted 
specimens were clamped at the center section and supported at the free 
end by an idler, in which the specimen was free to rotate. Also 
attached at the free end was the angle indicator and a moment arm. The 
moment arm was designed to act as a ten inch lever when loads were 
applied at a right angle to it through the eye bolt. A cable attached to 
the eye bolt fed through a dynamometer and pulley arrangement to a 
take-up ratchet. In this manner the applied moment was simply the 
dynamometer reading multiplied by the ten inch lever arm. 



Test II 
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As in Test I, an initial load was applied and all gages zeroed. 
Readings of strain, angular deformation, and applied torque were then 
taken in increments of applied torque. 



/ 
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4. 



Results and Discussion. 



The experimental value of Young's Modulus, as determined from 

Fig. 13 and Hooke's Law, was E = 10.1 x 10 psi. This is two percent 

lower than the handbook value'*' for the test material but was considered 

to be the most accurate figure available for that material sample. The 

2 

experimental E and the handbook values for the shear modulus, G, and 

3 

Poisson's ration, , were used as material constants in all 

calculations. An experimental value for the shear modulus could have 
been determined by a torsional test on the round stock prior to 
machining. From this test the shear modulus could be found from the 
torsional relationship for circular rods, 






T 

J Gr 



where J is the polar moment of inertia. The absence of such a test 
was an oversight in the test procedure. 

Table I lists the c>< and j3 coefficients for each of the 
three specimens, comparing the theoretical values with those determined 
by test. A graphical comparison of theoretical and experimental values 

^"United States Department of Defense, Strength of Metal Aircraft 
Elements: Military Handbook 5 (Washington: Government Printing 

Office, 1955) p. 114. 

2 Ibid. 

3 

Eshbach, 0. W. , Handbook of Engineering Fundamentals (New York: 
John Wiley & Sons, 1965) p. 12-56. 
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of the various coefficients relative to the initial twist rate is 
plotted in Fig. 14. It is significant to note the rather good 
degree of correlation attained by a relatively simple theory. The 

cross coupling coefficients, ©<■ ^ and ^21’ ft 12 an< * ft 21* 
were found to be equal as the analytical solution predicted. Prelim- 
inary investigation made it clear that membrane theory alone produced 
extremely erroneous values for °^22' m °difi ca tion introduced in 

the Analysis to obtain ^ f rom a theoretical ^ based on the 
theory of torsion proved to allow a reasonable prediction of this coef- 
ficient. It should be noted here that an error was introduced in the 
application of Chen Chu's modifying factor for the high initial rate of 
twist since he stated that there was an upper limit in his analysis, 
beyond which higher order effects were significant. £ 1 J This upper 
limit corresponds to the initial twist of Specimen #2, 



TABLE I 

TABULATION OF RESULTS 



Coefficient 


Specimen #1 
K = .2945 


Specimen #2 
K = .1473 


Specimen #3 
K = 0 


Theory 


Test 


Theory 


Test 


Theory 


Test 


11 


.866 


.892 


.963 


.947 


1.0 


.980 


12 


.0749 


.0668 


.0460 


.0389 


0 


0 


21 


.0749 


.0696 


.0460 


.0396 


0 


0 


22 


.0204 


.0169 


.0121 


.0115 


.00838 


.00897 


11 


1.69 


1.62 


1.27 


1.23 


1.0 


1.02 


12 


-6.21 


-6.40 


-4.82 


-4.16 


0 


.08 


21 


- 6.21 


-6.67 


-4.82 


-4.24 


0 


0 


22 


71.8 


85.5 


101.0 


101.3 


119.3 


111.4 
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The three experimental values of ^ ^ were obtained from Fig. 15, 
where e^ versus P* is plo tted for T* = 0 (Test I). The other coefficient 
obtainable from Test I was ^ the ratio of X to P* with T* =2 0. 

These values were obtained from Fig. 16. Likewise, Test II provided 
experimental determination for coefficients jS ^ an <3 [3 22 ’ 

(e /T*)p^_g and (X/T*)p^_Q, respectively, as obtained from Figs. 17 and 
18. These |3 coefficients were converted to experimental c< co- 
efficients by eqn. (37). 

Figs. 19, 20, and 13 are plots of strain data obtained from Test I 

for Specimens #1, #2, and #3 respectively. The strain reading from 

2 2 

every station, modified by the factor (1 + K j? ), was plotted versus 
2 

£ for P s 10,000 pounds, T s 0, This strain distribution across the 
chord is seen to plot in a linear fashion as predicted by eqn. (25), 
which is drawn as the solid line. Since the tangential strain at all 
chordwise stations was found to vary linearly with load, the value of 
10,000 pounds was simply a convenient basis and any other load would 
plot qualitatively the same. 

The strain distribution over the chord was essentially parabolic 
for twisted specimens, with the greatest strain realized at the mid- 
point of the chord. The strain relief evident at the outer fibers is 
directly attributable to the amount of angular deformation, or "un- 
winding" of these specimens. An interesting point to note in Fig. 19 
is that while the specimen is subjected to a 10,000 tensile load, the 
outer fibers are actually in compression. This is a result of the extreme 
angular deformation for the case of high initial twist rate. For the un- 
twisted specimen the strain distribution was uniform over the chord, as 
one would expect. 
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Figs. 21 through 23 are similar plots of strain distribution across 
the chord but for the condition T = 125 inch pounds, P * 0. Again the 
parabolic nature of the strain distribution is evident from the data 
for the twisted specimens, but there is no discernible pattern for the 
untwisted specimen. This last case appeared to be the only test in 
which the measured surface strain was not representative of the actual 
strain at the central fibers at each station. The "back to back" 
arrangement of strain gages on all specimens served as a check for this 
type of discrepancy. Test II of Specimen #3 was the only case in which 
a noticeable difference in strain readings was found from one side of 
the blade to the other. These differences were qualitative as well as 
quantitative, with e^ actually changing from tension to compression from 
one surface to the other. Even an average of each "back to back" pair 
of strain gages failed to provide a discernible pattern of behavior 
for this specimen. 

The reason for the lack of conformity of this test can most likely 
be traced to the test apparatus. The torsion idler ideally allowed the 
specimen to rotate completely free of friction but, at the same time, 
completely supported in all directions so no bending normal to the 
z-axis would occur. A compromise was adopted for this item in that the 
idler was simply a block support with a slightly oversize hole to allow 
free rotation. The resulting "play" in this hole allowed some degree of 
bending of all specimens as the moment was applied, but only in the case 
of Specimen #1, where the magnitudes of longitudinal strains due to 
torsion were so small, did the longitudinal strains due to bending 
become an important factor* 
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The untwisted specimen, of course, violates the assumption of 
“large" initial twist rates that was made in the analysis of membrane 
stresses. However, at low rates of twist the membrane stresses become 
insignificant in comparison with the shearing stresses in resisting an 
applied moment. The untwisted specimen was therefore included in the 
test as an important limiting case. It also served as a means of 
experimentally determining Young’s Modulus. 

The results of these tests have shown the theoretical approach 
herein utilized to be a good first approximation for defining the 
reaction of a twisted blade to a combined stress field. The results are 
not exhaustive, however, since specimens of many different initial twist 
rates and cross sectional aspect ratios would be required to establish a 
definitive pattern of response. An ability to apply completely arbi- 
trary combinations of load and moment in the test procedure would 
provide for a more thorough investigation of experimental-theoretical 
correlation. 

Second order effects which were neglected in this analysis assume 
much more importance at very low initial twist rates. It is assumed 
that membrane theory would prove more reliable for a thinner blade 
cross section; i.e., as the blade actually approached a pure membrane. 
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5 , 



Conclusions . 



The conclusions to be drawn from this paper are, for a twisted 
blade of narrow rectangular cross section: 

a. Modified membrane theory provides a good first approximation 
to the response of the blade to combined loads of tension and torsion. 

b. The strain distribution across the chord is essentially 
parabolic and is accurately predicted for the twisted specimens tested. 
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6 . 



Recommendations . 



It is recommended that further work in the area of twisted blades 
under combined loading be directed toward 

a. a more complete theoretical treatment, to include, for 
example, a detailed analysis of the shear stresses in the plane of 
the cross section as well as an investigation of the effect of higher 
order terms that were neglected in the present analysis, 

b. the development of an experimental apparatus that will 
allow a complete freedom of choice of combinations of the applied 
loads, and 

c. a confirmation of the validity of the conclusions made in 
this paper by tests on a number of specimens in which both cross 
sectional aspect ratios and initial twist rates are varied. 
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Fig. 1. Diagram of Analytical Model. 
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Fig. 2. The Blade Thickness, t. 



38 




39 



Strained Condition (line OP') 




Fig. 4. The Trs( Spocirrwns 
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Fig. 5. BLH Switching and Balancing Unit and Budd 
DATRAN Digital Strain Indicator 






Fi<'. 6. Display of Experimental Apparatus 
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4 3 



Experimental Set-up, Test 1 (view a) 




Fig, 8. Experimental Set-up, Test ] (view b) 
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Fig. 9. Experimental Set-up, Test I (view c) 
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Fig. 10. Experimental Set-up, Test II (view a) 
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Fig. 11. Experimental Set-up, Test II (view b) 
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Experimental Set-up, Test II (view c) 
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1000 




.8 .9 T.O 



Fig. 13. Strain Distribution Across the Chord, Specimen #3, 
I (P a 10,000 pounds). 
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Fig, 14. Comparison of Theoretical and Experimental Values of 
&<■ and A Coefficients Relative to Initial Twist Rate. 
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Fig. 15. Plot of ( e 0 / p *) x * 3 Q 
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Fig. 16. Plot of CX/P*) T * =0 
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avg 
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Fig. 18. Plot of (%/T*) p * =0 
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Fig. 19. Strain Distribution Across the Chord, Specimen #1, 
I (P = 10,000 pounds). 
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Fig. 20. Strain Distribution Across the Chord, Specimen #2, 
I (P a 10,000 pounds). 
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Fig. 21. Strain Distribution Across the Chord, Specimen #1, 
Test II (T * 120 inch pounds). 
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Fig. 22. Strain Distribution Across the Chord, Specimen #2, 
Test II (T a 129 inch pounds). 
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Fig. 23. Strain Distribution Across the Chord, Specimen #3, Test II. 
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APPENDIX A 



TABLE II 

TEST SPECIMEN SPECIFICATIONS 



o . 1 

Specimen 


Cross Section 


Blade 


d0/dz 


K 


0 „ 
max 




a 

(in. ) 


h 

(in. ) 


h 

avg 

(in.) 


Length 
(in. ) 


(rad/in. ) 




#1 


.75 


.200 


.182 


9.3 


.3927 


.2945 


.2864 


#2 


.75 


.200 


.196 


9.2 


.1963 


.1473 


.1473 


#3 


.75 


.200 


.200 


12.7 


0 


0 


0 



* T = R/h = 3.125 for all specimens 

£ 

\ 
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TABLE III 



STRAIN GAGE PLACEMENT 
SPECIMEN #1 



Gage 1 

Number 


2 

Gage 

Type 


„ 3 

Gage 

Factor 


<£ 

Distance 


f‘ 


1 


A 


1. 


99 


0.563 


0.563 


2 










.563 


.563 


3 










0.000 


0.000 


4 










0.000 


0.000 


5 










.594 


.626 


6 










.563 


.563 


7 










.563 


.563 


8 










.563 


.563 


9 










0.000 


0.000 


10 










0.000 


0.000 


11 










.563 


.563 


12 










.594 


.626 


13 










.344 


.210 


14 










.344 


.210 


19 










.375 


.250 


20 










.375 


.250 


25 










.344 


.210 


26* 










.375 


.250 


31 










.375 


.250 


32 


1 


f 


> 


r 


.375 


.250 


37 


B 


2. 


10 


.750 


1.000 


38 










.750 


1.000 


39 










.750 


1.000 


40 




f 


1 


r 


.750 


1.000 



Listed so the "back to back" gages appear as a pair. 

"A" corresponds to type A-7. "B" corresponds to type 

FAP 12-12 (S13). 

All gage resistances are 120 - .5 ohms. 

Faulty gage. 



61 



TABLE IV 



STRAIN GAGE PLACEMENT 
SPECIMEN #2 



Gage 1 

Number 


2 

Gage 

Type 


Gage 3 

Factor 


£ 

Distance 




1 


A 


1.99 


0.58 


0.598 


2 


A 


1.99 


.54 


.518 


3 


B 


2.09 


0.00 


0.000 


4 


B 


2.09 


0.00 


0.000 


5 


A 


1.99 


.55 


.538 


6 


B 


2.09 


.56 


.558 


7 


B 


2.09 


.38 


.257 


8 


B 


2.09 


.35 


.218 


9 


A 


1.99 


.36 


.230 


10 


B 


2.09 


.37 


.243 


11 


A 


1.99 


.34 


.205 


12 


B 


2.09 


.40 


.284 


13 


B 


2.09 


.39 


.270 


14 


B 


2.09 


.35 


.218 


15 


A 


1.99 


.55 


.538 


16 


A 


1.99 


.59 


.619 


17 


B 


2.09 


0.00 


0.000 


18 


B 


2.09 


0.00 


0.000 


19 


A 


1.99 


.55 


. .538 


20 


A 


1.99 


.55 


.538 


21 


B 


2.10 


.75 


1.000 


22 


B 


2.10 


.75 


1.000 


23 


B 


2.10 


.75 


1.000 


24 


B 


2.10 


.75 


1.000 



''"Listed so the "back to back" gages appear as a pair. 

2 

"A" corresponds to type A-7. "B" corresponds to type 

FAP 12-12 (S13 ) . 

All gage resistances are 120 - .5 ohms. 
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TABLE V 



STRAIN GAGE PLACEMENT 
SPECIMEN # 3 



Gage 1 

Number 


Gage 

Type 


2 

Gage 

Factor 


Distance 




3 


FAP 12-12 


2. 


09 


0.00 


0.00 


9 










0.00 


0.00 


2 










.35 


.22 


10 










.36 


.23 


4 










.36 


.23 


8 










.35 


.22 


1 










.56 


.56 


11 










.55 


.54 


5 










.57 


.58 


7 










.56 


.56 


6 3 










.75 


1.00 


12 


> 




'I 


f 


.75 


1.00 



1 

Listed so the "back to back" gages appear as a pair. 

2 

All gage resistances are 120 - .5 ohms. 

3 

Faulty gage. 
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TABLE VI 



TEST I DATA, SPECIMEN #1 



P (pounds) 

P* x 10 3 
Y ' ( 1/16 inch) 1 
X x io 3 


2000 

.66 

10 

4.4 


4000 

1.32 

20 

8.8 


6000 

1.98 

30 

13.2 


8000 

2.64 

40 

17.6 


10000 

3.30 

50 

22.0 


Gage Number 




S train (microinches 


per inch) 




i 


350 


715 


1080 


1460 


1860 


2 


325 


670 


1020 


1400 


1780 


3 


1100 


2200 


3280 


4370 


5440 


4 


1075 


2150 


3220 


4280 


5340 


5 


280 


580 


880 


1200 


1530 


6 


350 


720 


1100 


1490 


1890 


7 


325 


650 


1000 


1350 


1720 


8 


325 


660 


1010 


1370 


1750 


9 


1020 


2030 


3040 


4040 


5040 


10 ' • 


1120 


2230 


3350 


4460 


55 70 


11 


300 


600 


920 


1250 


1600 


12 


310 


630 


960 


1310 


1670 


13 


860 


1730 


2600 


3490 


4380 


14 


830 


1650 


2480 


3320 


4160 


19 


810 


1620 


2430 


3260 


4090 


20 


750 


1520 


2290 


3080 


3880 


25 

26 * 


820 


1640 


2490 


3350 


4210 


31 


740 


1490 


2250 


3040 


3820 


32 


800 


1620 


2440 


3270 


4100 


37 


-190 


-350 


-490 


-600 


-700 


38 


-150 


-280 


-390 


-480 


-540 


39 


-170 


-320 


-460 


-570 


-660 


40 


-200 


-380 


-540 


-670 


-790 


e x 10 3 
o avg 


1.08 


2.15 


3.22 


4.29 


5.35 



^See eqn. (48). 
2 

Faulty gage. 
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TABLE VII 



TEST I DATA, SPECIMEN #2 



P (pounds) 


2000 


4000 


6000 


8000 


10000 


P* X 10 3 


0.66 


1.32 


1.98 


2.64 


3.30 


Y' (1/16 inch) 1 


6.40 


12.8 


19.2 


25.2 


31.2 


X * 10 3 


2.85 


5.70 


8.54 


11.21 


13.88 


Gage Number 




Strain (microinches 


per inch) 




i 


520 


1060 


1610 


2160 


2720 


2 


560 


1120 


1680 


2250 


2820 


3 


810 


1620 


2440 


3245 


4060 


4 


820 


1630 


2450 


3265 


4080 


5 


510 


1040 


1570 


2130 


2680 


6 


535 


1080 


1650 


2240 


2840 


7 


685 


1370 


2070 


2 760 


3460 


8 


700 


1400 


2080 


2 790 


3470 


9 


650 


1310 


1980 


2660 


3345 


10 


700 


1400 


2120 


2830 


3560 


11 


700 


1400 


2100 


2800 


3480 


12 


700 


1400 


2120 


2830 


3550 


13 


770 


1535 


2310 


3070 


3820 


14 


750 


1490 


2240 


3000 


3740 


15 


550 


1110 


1680 


2250 


2830 


16 


490 


1000 


1520 


2040 


2570 


17 


840 


1680 


2530 


3370 


4220 


18 


770 


1560 


2350 


3140 


3930 


19 


570 


1160 


1740 


2330 


2910 


20 


530 


1075 


1620 


2190 


2 760 


21 


400 


800 


1235 


1670 


2120 


22 


350 


720 


1120 


. 1540 


1970 


23 


340 


710 


1100 


1510 


1930 


24 


385 


750 


1190 


1620 


2055 


e x 10 3 
o avg 


.81 

1 


1.62 


2.44 


3.26 


4.07 





1 



S 



ee eqn. 



(48). 
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TABLE VIII 



TEST I DATA, SPECIMEN #3 



P (pounds) 

P* x lo 3 
^ (1/16 inch)* - 
X x 10 3 


2000 

.66 

0 

0 


4000 

1.32 

0 

0 


6000 

1.98 

0 

0 


8000 

2.64 

0 

0 


10000 

3.30 

0 

0 


Gage Number 




S train (microinches 


per inch) 




3 


650 


1320 


1990 


2670 


3340 


9 


675 


1355 


2030 


2710 


3400 


2 


650 


1310 


1975 


2650 


3320 


10 


665 


1355 


2030 


2710 


3390 


4 


650 


1305 


1970 


2650 


3330 


8 


660 


1330 


2000 


2680 


3360 


1 


650 


1320 


2000 


2690 


3370 


11 


665 


1345 


2020 


2 700 


3380 


5 


645 


1295 


1955 


2630 


3310 


7 

6 2 


665 


1345 


2010 


2 700 


3380 


12 


665 


1345 


2015 


2700 


3380 


e x 10 3 
o avg 


.66 


1.34 


2.10 


2.69 


3.37 



^See eqn. (48). 
2 

Faulty gage. 
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TABLE IX 



TEST II DATA, SPECIMEN #1 





T (inch pounds) 


30 


46 


85 


120 


150 


180 


210 


T* x 10 6 


13.1 


24.4 


37.0 


52.2 


65.3 


78.3 


91.4 


y' (1/16 inch)* 


2.7 


4.9 


7.3 


10.1 


12.8 


15.1 


17.7 


X x 10 3 


1.19 


2.17 


3.21 


4.44 


5.63 


6.64 


7.79 


Gage Number 




Strain (microinches per 


inch) 




i 


100 


181 


273 


373 


463 


545 


635 


2 


98 


180 


267 


370 


460 


540 


625 


3 


-95 


-170 


-256 


-349 


-430 


-510 


-590 


4 


-88 


-158 


-242 


-332 


-410 


-485 


-565 


5 


105 


198 


303 


419 


520 


615 


715 


6 


83 


158 


240 


329 


410 


485 


560 


13 


-21 


-41 


-70 


-95 


-123 


-144 


-169 


14 


-33 


-62 


-86 


-123 


-154 


-183 


-212 


19 


-10 


-16 


-26 


-34 


-42 


-50 


-56 


20 


-18 


-25 


-33 


-45 


-56 


-67 


-80 


37 


221 


409 


620 


855 


1060 


1255 


1455 


38 


214 


398 


600 


825 


1020 


1210 


1400 


e x 10^ 
o avg 


-92 


-164 


-249 


-341 


-420 


-497 


-577 



1 



S 



ee eqn. 



(48). 
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TABLE X 



TEST II DATA, SPECIMEN #2 



T (inch pounds) 
T* x 10 6 

(1/16 inch)l 
X x 10 3 


24 

10.56 

2.4 

1.07 


49 

21.56 

5.3 

2.36 


75 

33.00 

7.8 

3.47 


104 
45.76 
10.6 
4. 72 


129 

56.76 

13.3 

5.92 


153 

67.32 

15.3 

6.81 


182 

80.08 

17.9 

7.97 


Gage Number 




Strain (microinches per 


inch) 


" 

i 


ii 


-15 


-34 


-52 


-74 


-87 


-104 


-123 


12 


-1 


-8 


-13 


-18 


-21 


-26 


-2 7 


- i 13 


-25 


-48 


-70 


-97 


-121 


-142 


-167 


14 


-15 


-28 


-45 


-65 


-76 


-91 


-104 


15 


41 


76 


113 


152 


186 


216 


247 


16 


39 


82 


121 


163 


202 


237 


2 72 


17 


-58 


-112 


-168 


-231 


-283 


-332 


-387 


18 


-41 


-79 


-118 


-161 


-198 


-230 


-268 


19 


28 


54 


81 


108 


132 


153 


176 


20 


41 


89 


127 


179 


220 


257 


300 


23 


115 


229 


343 


468 


575 


670 


778 


24 


112 


223 


333 


457 


559 


652 


755 


e x 10 
o avg 


-50 


-96 


-143 


-196 


-240 


-281 


-328 



1 



See eqn. 



(48). 
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TABLE XI 



TEST II DATA, SPECIMEN #3 



T (inch pounds) 
T* x 10 6 
y (1/16 inch) 1 

X x io 3 


16 

7.0 

2.4 

.77 


33 

14.5 

5.2 

1.67 


49 

21.6 

7.5 

2.42 


68 

29.9 

10.4 

3.35 


86 

37.8 

13.4 

4.31 


107 

47.1 

16.3 

5.25 


126 

55.4 

19.1 

6.15 


Gage Number 




Strain (microinches per inch) 




3 


-2 


-6 


-15 


-21 


-25 


-30 


-36 


9 


5 


10 


16 


23 


29 


36 


45 


2 


-1 


-4 


-7 


-15 


-15 


-20 


-24 


10 


12 


27 


41 


57 


67 


84 


99 


4 


-12 


-23 


-38 


-50 


-64 


-79 


-93 


8 


-1 


-2 


-6 


-8 


-12 


-11 


-12 


1 


-4 


-10 


-14 


-23 


-31 


-38 


-44 


11 


17 


30 


48 


63 


77 


93 


110 


5 


-21 


-40 


-68 


-92 


-117 


-147 


-171 


7 

6 2 


-5 


-10 


-14 


-18 


-21 


-26 


-31 


12 


12 


18 


27 


33 


38 


43 


48 


e x 10^ 
o avg 


1.5 


2.0 


0.5 


1.0 


2.0 


3.0 


4.0 



^See eqn. (48). 
2 

Faulty gage. 
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APPENDIX B 



Calculation of Young's Modulus, E. 

From Fig. 13: 

P = 10,000 pounds 

= 336 x 10 inches/inch 

avg 



E 




P 

2ahe t 



10000 

2(.75)(.2)(. 000336) 



10.1 



* io 6 



psi 



Calculation of o< Coefficients, 

The calculation of the first three c* coefficients involves the 
solution of the appropriate integrals in eqns . (29) and (30). These 
can be broken down into basic integrals of the form 





(l + K x f l ) L 




where n = 0, 2, and 4. These integrals will be referred to as 1^, I 
and I^, respectively. 



I 







df 



I 

2.(1 + K 2 ) 



+ 



I 

Tk 



fan ' 1 K 



but 



1 + K‘ 



= 1 



K 2 + 



K 



K 



+ 



and tan K = K 





5 



+ 



Substituting, 



x i = icV- 0] 

x i = l - j K a +• f- K 4 + "* 

f' ? l - 1 i 

11 'Jo U + tff*) 3 - ^ " II^TTTh^T + if 5 * ar * ' K 

1 2 “ t(/ "IT 1 + 1 ’ K' i + K 4 ^'*) +■ (-fat - -j + JH’ 2 - 

X 2= J"# ... 

r' i j 

13 = J 0 l n^yY = ifViK*) ' "h 7 2 * 

x 3 - -^[o-* 2 ** 4 -***-) -(/- j **+4 * 4 - ^ 






1 






From eqns. (29) and (30), 



°<n = I, - r = (1 K' 4 * 4 -) k 4 +- 

°<11 = I - ^'^ + t')K' i 4j(jf2r)K 4 f <** 

®<12 * **, s 

*12 = (3 K*V~) k* r +") 
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>11- 



°<12 “ 3"K - ^ (2 + r) K J f (3 4-zr) - 



From eqn. (47), 



°^22 



£ h«t<j (.C>6>7g. -.21 h ^)[/ +/%( I +/*■)( &*naie T^Z ) ] + 

ZcJE 



2 

(47) 

°<#i 



Substituting into the above three expressions for c-C^, anc * 

o< 22 , there is 



°Vi 55 



°*11*2 



°<11>3 = 



°^12 ) 1 * 



°^12' ) 2 ” 



1 - .148 + .014 * .866 



1 - .037 = .963 



1 



.0982 - .0262 + .0029 = .0749 



.0491 - .0033 = .046 



12 >3 



= 0 



°^22 ) 1 



5,9 X/o° (./ ?2.) 
Z (. 70 5 10.1 K/0 






[l <■&(!*■> S)(. t K 4 j£) 2 ]J 



Zl (, HZ )][ 

(.0741) *" 

, UeL 
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o^ 2 2 = .01516( .5 - , 0382)(1 + .988) + .0065 
* * 0139 + '° 065 85 * 0204 
o<' 22 ) 2 = - 0176 (- 5 ' .0412 )( 1 + .222) + .0022 
<* 22 ) 2 = .00986 + .0022 = .0121 
0 ^ 22)3 = .0183( .5 - .042) (1) + 0 = .00838 



Calculation of coefficients. 

The expressions for fB coefficients in terms of « 
ficients are to be found in eqns. (33) and (34). Let 



D = ex' 11 ~ 



(2 



11 ' 



P 12 



0 

" 2 






22 



D x = ( .866 )( . 0204 ) - (.0749) = .0121 



P 11*1 = ( .0204)/( . 0121 ) = 1.69 



coef- 
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ft I2h = ( -.0749 )/ ( . 0121 ) = -6.21 



^ 22 } 1 = ( .866)/ ( .0121) = 71.8 
D 2 = ( . 963 ) ( .0121) - ( .046) 2 = .00954 
|§ n ) 2 = (v0121)/(. 00954) = 1.27 
|S i 2 ) 2 « (-.046)/(. 00954) = -4.82 
(£ 22 ) 2 = (.963)/(. 00954) = 101 
D 3 = (1) (.00838) = .00838 
^ n ) 3 = (.00838)/(. 00838) = 1.0 

ft 12 ^3 = ° 

(S 22 ) 3 = (.963)/(. 00838) = 119.3 



Calculation of theoretical strain distribution. 



The theoretical strain distribution across the chord is given 
by eqn. (25), 



e t (i + k ¥> - % + * k / 2 



(25) 
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a) Test I: P = 10,000 pounds, T = 0 

From eqns. (35) and (36) 

e o - &U p * 

X = 0 21 ** 

Specimen #1: 

e Q = 1.69(3300 x lO -6 ) = 5577 x lo" 6 

X = -6.21(3300 x 10" 6 ) = -20500 x 10" 6 

e t (l + K 2 ^ 2 ) * 5577 x 10" 6 - 20500 x 10 _6 ( .2945) j* 1 

e t (l + K 2 ^ 2 ) = 5577 x 10' 6 - 6035 x lO -6 ^^ 

Specimen #2: 

e Q = 1.27(3300 x 10 _6 ) = 4191 x lO* 6 

X = -4.82(3300 x 10" 6 ) = -15,906 x 10‘ 6 

e t (l + K^ 2 ) = 4191 x 10 -6 - 15,906 x 10~ 6 ( . 1473) ^ 

e c (l + K 2 £ 2 ) = 4191 x 10“ 6 - 2343 x 10 _6 jf 2 
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Specimen #3: 



e Q = 1.0(3300 x 10" 6 ) = 3300 x 10 _6 

X - 0 

e t (l + K 2 ^ 2 ) = 3300 x 10” 6 

b) Test II: T ^ 125 inch pounds, P = 0 

From eqns . (35) and (36) 

% * ( 3 12 T * 

X - (3 22 T * 

Specimen #1: 

e Q = -6.21(52.2 x 10” 6 ) = -324 x 10 -6 

X = 71.8(52.2 x 10" 6 ) . 3748 x 10 -6 

e t (l + K 2 ^ 2 ) = -324 x 10 -6 + 3748 x 10 _6 ( .2945 ) f 2 

e t (l + K 2 ^ 2 ) = -324 x 10 -6 + 1103 x 10~ 6 J 2 
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Specimen #2: 



e = -4.82(56.76 x 10" 6 ) = -273 x 10‘ b 
o 

%, = 101(56.76 X lo" 6 ) = 5733 x 10* 6 

e t (l + K 2 ^ 2 ) = -273 x 10' 6 + 5733 x 10 

e t (l + K 2 £ 2 ) = -273 x 10' 6 + 844 x 10 _ ' 

Specimen #3: 

e = 0 
o 

X = 119.3(55.4 x 10' 6 ) = 6610 x 10' 6 
e t (l + K 2 ^ 2 ) = 0(6610 x 10" 6 )^ 2 = 0 



(.1473)^ 2 

,2 
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